Spin wave excitations in the antiferromagnetic 

Heisenberg-Kondo model for heavy fermions. 

oo ■ 

o: 

^ ■ M. Acquarone^ 

3 ; 

^ . C. I. Ventura^ 

^^ , ^IMEM-CNR and Dipartimcnto di Fisica, Universita di Parma, 43100 Parma, Italy, 

r— 1| ^Centro Atomico Bariloche, 8400 - Bariloche, Argentina , 

^ , E-mail: acquaroneOf is .unipr . it, venturaScab. cnea.gov. ar 

U: 

. ' Abstract. 

C^ I Recent inelastic neutron scattering experiments in Celna and CePd2Si2 single 

C I crystals measured spin wave excitations at low temperatures. These two heavy fermion 

compounds exhibit antiferromagnetic long-range order, but a strong competition 
between the Ruderman-Kittel-Kasuya-Yosida(RKKY) interaction and Kondo effect is 
qI ' evidenced by their nearly equal Neel and Kondo temperatures. Our aim is to show how 

O/ magnons such as measured in the antiferromagnetic phase of these Ce compounds, can 

be described with a microscopic Heisenberg-Kondo model introduced by J.R.Iglesias, 
C.Lacroix and B.Coqblin, used before for studies of the non-magnetic phase. The 
model includes the correlated Ce-4/ electrons hybridized with the conduction band, 
^h I where we also allow for correlations, and we consider competing RKKY (Heisenberg- 

O^l ' like Jh) and Kondo (Jk) antiferromagnetic couplings. Carrying on a series of 

" . , unitary transformations, we perturbatively derive a second-order effective Hamiltonian 

C^^ I which, projected onto the antiferromagnetic electron ground state, describes the spin 

^^ • wave excitations, renormalized by their interaction with correlated itinerant electrons. 

f^ , We numerically study how the different parameters of the model influence the 

ILJ ' renormalization of the magnons, yielding useful information for the analysis of inelastic 

neutron scattering experiments in antiferromagnetic heavy fermion compounds. We 
also compare our results with the available experimental data, finding good agreement 
j^ ■ with the spin wave measurements in cubic Celus. 
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1. Introduction 

The description of heavy fermion compounds is challenging due to the rich variety 
of phase diagrams they present, and the anomalous physical properties which may 
be found. Among them, appear the Ce and U compounds which exhibit long-range 
antiferromagnetism (AF) at low temperatures (for example, among antiferromagnetic 
Ce compounds, CeRh2Si2 exhibits the highest ordering temperature T]^] = 36K, 
with local magnetic ordered moments of 1.34 — 1.42/^^ per Ce, i.e. relatively large, 
compared to the full Ce^+ free-ion value: 2.54 /iB).p!] Depending on the particular 
compound, [2] the antiferromagnetism takes different forms (magnitude of the local 
moments varies widely: e.g. 0.001 fiB as in CeRu2Si2 or 0.02 /i^ in UPts, to 1.55 
fiB as in UCus; as does the spin configuration: with three-, two- or one- dimensional 
AF structures observed). Antiferromagnetism may also appear competing or coexisting 
with superconductivity, [31 H] for which spin fluctuation-mediated pairing mechanisms 
are explored. [5] Non- Fermi liquid behaviour may appear, [6] and quantum criticality 
has become a subject of intensive study in these compounds, both experimentally and 
theoretically. [H [TJ [H [9l [TOl [11] The crossover from the antiferromagnetic state to the 
non-magnetic heavy fermion state, which can be tuned by pressure, doping or magnetic 
field, is one of the most interesting problems in strongly correlated /—compounds. 

The physical properties of these compounds are determined by the strongly 
correlated /—electrons present (4/ in Ce; 5/ in U) and their hybridization with 
the conduction band. The RKKY indirect exchange interaction between /— local 
magnetic moments, favouring the establishment of long-range magnetic order, competes 
with the screening of these moments by the conduction electrons, described by 
the Kondo effect. [12] This competition is the subject of the Doniach diagram. [T3] 
which compares the variation of the Neel and Kondo-impurity temperatures with 
increasing antiferromagnetic intrasite exchange coupling Jk, between local /-moments 
and conduction electron spins. Compounds with similar magnetic ordering temperature 
Tjv and Kondo temperature T^, the temperature below which magnetic susceptibility 
saturates indicating coherent Kondo-singlet formation, are ideally suited to the study 
of this RKKY-Kondo competition. In this regard, experiments on CeM2Sn2 (M=Ni, 
Ir, Cu, Rh, Pd, and Pt: Ad or 5d transition metals) [H] and CeX2Si2 (X=Au, Pd, Rh, 
Ru) [15] were undertaken, indicating that departures between theory and experiments 
resulted from the use of Kondo impurity relations. The Kondo-lattice model, instead, 
consisting of a lattice of local magnetic moments coexisting with a conduction band, 
has proved appropriate for the description of many 4/ and 5/ materials, in particular 
most Ce (or Yb) compounds, respectively corresponding to a configuration close 
to 4/^ (or 4/^^), where one 4/ electron (or hole) interacts with the conduction 
electrons. [TOl [T6| \T7\ [T8] In 1997 a revisited Doniach diagram was introduced, including 
short-range antiferromagnetic correlations in the Kondo lattice, in order to improve the 
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description and, in particular, to account for the observed pressure dependence of Tk 
in CeRh2Si2.[16] The situation is more complex in Uranium compounds, where U has 
a 5/" configuration with n=2 or 3, since the 5/ electrons are much less localized than 
the 4/ electrons of rare earths. Regarding spin dynamics, it is not clear that Ce and U 
compounds are intrinsically similar, p] In the following we will focus on Ce-compounds, 
except otherwise specifically stated. 

Experimentally, while the magnetic response due to Kondo spin fluctuations in the 
paramagnetic state of heavy fermions is well studied, relatively little is known about 
the nature of the magnetic excitations in the ordered phase of Kondo lattices, [2] on 
which our present study will focus. Being still an unsolved problem how to describe on 
equal terms both the Kondo effect and antiferromagnetism. |19j here we will focus on 
systems with relatively large local magnetic moments and study them deep inside the 
antiferromagnetic phase: far from the antiferromagnetic quantum critical point, where 
spin fluctuations would become more relevant. 

A few years ago CePd2Si2[2n] single crystals were studied with inelastic neutron 
scattering: below the antiferromagnetic ordering temperature strongly dispersive spin 
wave excitations were found, with an anisotropic damping, which coexisted with the 
Kondo-type spin fluctuations also present above T/y. At T = 1.5K these spin waves were 
measured along various BZ paths: they were found to present an energy gap of 0.83 meV 
and to extend up to almost 3.5 meV. CePd2Si2 has a bcc tetragonal structure, and its 
antiferromagnetic ground state is characterized by propagation vector: q= (1/2, 1/2, 0), 
with ordered moments: S = 0.66/iB, Tn = 8.5K and Tk = lOK, and linear electronic 
specific heat coefficient 7 = 250 mJ/mol K^. Under pressure application, at 28.6 
kbar the system undergoes a transition into a superconducting phase with critical 
temperature of 430 mK.[2T| [22] More recently, inelastic neutron studies of Celus single 
crystals were performed, with similar results. [21] Well defined spin wave excitations 
with a bandwidth of 2 meV and a gap of 1.28 meV were found in the antiferromagnetic 
phase, |2T] coexisting with Kondo-type spin fiuctuations and crystal- field excitations 
which also appeared above T/v = 10K= Tk- Celns crystallizes in a cubic (fee) 
structure, with an antiferromagnetic structure characterized by magnetic propagation 
vector q = (1/2,1/2,1/2), with ordered moments: S = 0.5fj,B and 7 = 130 mJ/mol 
K^. Under application of pressure, at 26.5 kbar the system undergoes a transition 
into a superconducting phase with critical temperature of 200 mK.[2T] Recently, similar 
antiferromagnetic magnon excitations were measured by inelastic neutron scattering in 
CeCu2, [23] an anisotropic antiferromagnetic heavy fermion compound with Tjy = 3.5K 
and Tk = 4K. 

In next section, we will briefly introduce the microscopic Heisenberg-Kondo model 
proposed by J.R.Iglesias, C.Lacroix and B.Coqblin,[T6] to study the non-magnetic phase 
of heavy fermion AF compounds, to which we shall add conduction electron correlations. 
We will then present our calculation for the renormalization of spin wave excitations due 
to their interaction with the correlated conduction electrons (the Appendix complements 
this section). In Section^ we will discuss the results of our study, show how the different 
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parameters of the model influence the renormalization of the magnons, and compare 
our results with the available experimental data. [201 121] In Section [6] we summarize and 
point out that the present work should yield useful information for the analysis and 
prediction of inelastic neutron scattering experiments in heavy fermion AF compounds, 
as CeRh2Si2. 

2. Microscopic model, and perturbative approach. 

In order to describe the Ce-heavy-fermion systems exhibiting antiferromagnetic long- 
range order we have used the microscopic model which has been proposed by Coqblin et 
al.[16j to describe the competition between the Kondo effect and the RKKY interaction, 
in compounds where departures from the original Doniach picture |13] appear. In 
principle, both the RKKY magnetic coupling and the Kondo effect can be obtained 
from the Kondo intrasite-exchange term, but when dealing with approximations it 
is difficult to assure that both effects are taken into account if an explicit intersite 
exchange (as the effective RKKY interaction or, depending on the system, also the direct 
exchange) is not included in the Hamiltonian. [TO] The mo del [16] consist s of a Kondo 
lattice, featuring local magnetic moments coupled both to conduction electrons, by a 
Kondo-type interaction J^ , and, among themselves, by an antiferromagnetic RKKY- 
type exchange J^j > 0. The moments are assumed to order below the Neel temperature 
T/v, but we will concentrate on the zero-temperature limit. The Hubbard-correlated 
conduction electrons occupy a non-degenerate band. Therefore, the model may be 
represented in standard notation by the following Heisenberg-Kondo Hamiltonian: 

H = Hueis + H^^^^ + H Kondo (l) 

Hueis = Jr 2_^ ^i'^m — Ba I ^ S^ — ^ S^j 1 (2) 

i{m) KlaA jeB J 

Hband = Yl ^^^*- + Yl ^im4aCma + ^^ J] ^^^^H (^) 

iff i{m)a i 

HKondo = Jx I $^S, ■ S, + 5^S^. • S^. J (4) 

\l&A jeB J 

where we have included in Hheis an anisotropy field Ba > 0. which is physically realized 
by the crystal field in the heavy fermion systems. It will be shown to play a crucial role 
with respect to the stability of the AF spin waves. 

Here i (m) indicates that the lattice site index m runs over the z nearest neighbours 
of site i which, in turn, runs over the A^ sites of the full lattice. The itinerant electron 
spin is S;/-), while A,B label the two interpenetrating sublattices with A^/2 sites each, 
where local moments, Si^_^ and 5*^^^ , with opposite moment direction, sit. The Kondo 
exchange Jk could in principle have either sign (though for heavy fermion compounds, 
it would be antiferromagnetic: J^- > 0). 
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2.1. Diagonalization of Hfjeis- 

We diagonalize the Heisenberg term, by representing the local moments operators in the 
Holstein- Primakoff approximation, which is appropriate at temperatures much lower 
than T^r . Namely, we take 



Si — [Si , Si ,Si) 

2Sbi,V2SblS-blbi) 

2Sb,:,-iS-b\bj) 



5,~ 



2Sb], 



(5) 



'(t) 



.(t) 



where b} and b, are the bosonic spin-deviation operators in sublattices A and B, 
respectively. Going to reciprocal space in the reduced Brillouin zone (RBZ) one gets 
an anharmonic Hamiltonian in the spin-deviation |6j, 6^} operators. It's convenient for 
what follows to recall that a generic anharmonic Hamiltonian of the form: 

<? 1 

where Fg, Gq are c-numbers, is diagonalized by a Bogolyubov transformation, which 
introduces the AF spin wave operators {^a^^gja^^ according to: 



i^an/^ = 5^i^,(6j^ + 6L,6_, 



bl 



JCh (^J + a.^Sh (^J 
^ Sh (^J + a^Ch (^J 



a 



The diagonalization condition is Th (2'i9 ) = —Gq/Fq 
Hamiltonian (defining sgn {x) = x/ \x\ ) reads: 



(7) 
The diagonalized 



e^HanhG ^ = Hdiag = ^ sgu (Fg) ^ F^ - G^ (^al^a^^ + aja^ j 



(8) 



Stability of the system requires the renormalized frequency to be real and positive. 
Reality imposes the condition F^ > G^, i.e. the anharmonic part must have an 
amplitude smaller than the harmonic part. If the renormalized frequencies are real, 
their positiveness is assured by the additional constraint Fq > 0. 

In our case the diagonalization condition for Hhcis reads: 



Th (2^g 



Gn 



JjjS 



Z U TT O 



B„ 



E cos {q\ 



(9) 



In the absence of the interaction with the fermions, i.e. in the limit of vanishing Jk, 
the frequency of the bare AF spin waves would be [z is the number of nearest neighbors 
and A;j is the vector joining two n.n. sites) 



H 



Heis 



hn„ 



E^^. 



S«g + 2 



ZjrrS + B, 



a) \ 



'^H'5EA„cos(gAij) 



-\ 2 



zJj^S 



B„ 



(10) 



;ii) 
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2.2. Diagonalization of Hf,and- 

The band electrons described by Hband can be in either the paramagnetic (PM) or AF 
state, according to the values of the bandwidth W = 2zt and of the Hubbard correlation 
U . To diagonalize Hf,and we will use a reformulation of Gutzwiller's variational approach 
for the description of antiferromagnetism in narrow bands due to Spalek et al. [21] This 
approach allows to connect smoothly the PM state for U <^ W to standard (mean- 
field) Slater band-insulator ioT U ^ W and to the localized Mott antiferromagnetic 
insulator for U S> W. One expresses the correlation-induced bandwidth reduction in 
the paramagnetic (PM) state by a Gutzwiller-type factor ^{n,p) depending on the 
band filling n and on the probability of double occupancy p = N~'^ J2i {^ii^ii) ■ The 
correlated band energies e^ are given by: 

s'',=^n,p)el (12) 



^{n,p) = l , . . ^ 2 

2 — nj \ n^ 

where e^ are the uncorrelated band energies. The U- depending optimal value of p is 
found at zero temperature by minimizing the PM energy 

E^^^ = Hn,p)J2el{n,^) + NUp (13) 

ka 

at given n and U . 

Assuming that the electrons have an AF ground state of Neel-type, one adopts the 
standard Slater formalism only with the PM energies renormalized according to EqIT2l 
We therefore introduce the fermion operators for this AF Slater-type state < c^ka^Pka \ 
by the transformation 

cL = "L cos Cka - PL sin Cka 

4+s,a = "L sin Cka + PL cos Cka (14) 

where Q =(1/2, 1/2, 1/2) in units of 27r/a is the wavevector characterizing the AF 
magnetic state and fc is a wavevector belonging to the reduced Brillouin zone (RBZ) 
defined by |k| < \Q\. The diagonalization condition for Hband in a lattice with an 
inversion center yields: 

tan (2^.) = -^^7^ = 0- tan (2Cfc) , (15) 



[s) 



afS 



of the itinerant electrons 



where the amplitude of the AF order parameter 

(staggered magnetization, or band electron polarization), is given by: 

(^)./3| = ^ E (K.)-(^L))sin2a (16) 



2A^ 

keRBZ,! 



Let's stress that EqUnidoes not set the sign of (s)^^. When introducing the interaction 
with the ordered local moments, its sign will be set according to 

(s)a/3 = -Sgn ( Ji^) (s)„/3 (17) 
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so that a positive Jk means an antiparallel orientation with respect to the local moment 
staggered magnetization, assumed positive by definition, and viceversa. 
The resulting diagonal bare electron Hamiltonian is: 



H, 



AF 
band 



E 



Etnt + Elnl 



UN 



-{s)l 



keRBZ,(T 

where the bare electron eigenenergies (actually spin- independent) read (x = a,/3): 

EL = \un + (1 - 25.„) ^$%2 + f/2 (,)2^, (19) 

where 5xa denotes the Kronecker delta. Notice that a is the lower subband. The 
subband filling factors {nL) = {^l -a) — l^^P (E'L — fi) /ksT + 1]~ also depend on 
(s)^^ through E^^ so that Eq{T6]has to be solved self-consistently. 



3. The Kondo coupling. 

Taking into account the orientation of the local-moment magnetization, we shall 
distinguish in the Kondo term the longitudinal H^ from the transverse H^ 
contributions. By taking as positive z direction the direction of S^^^^ , they are defined 
as: 






^ 2 



.leA j&B 



(20) 
(21) 



3.1. The longitudinal part H^ of the Kondo coupling term. 



It is convenient to rewrite the longitudinal Kondo term by decomposing it into the two 
sublattices contributions. For a given sublattice X = A,B we have from Eq[ 






IGX 



= H^^, + H%^ {X = A, B) (22) 

This expression, rewritten in terms of the operators < Op > and < apa,l3pj \ which 
diagonalize, respectively, Hueis and Hi^and , yields several contributions. The first one, 
not containing Bose operators, is: 



Ja 



H'M + Hk = -^SJ2[{ 



"^pa "^pa 



)sin(2Cp)] 



p.a 



J. 



+ ~V^ E ^ ^°^ (^^P) {Pla(^p<r + (^Ijpa) 



(23) 



p,a 
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Notice that the second term in Eq|23] describes a Kondo-induced hybridization between 
the bare electrons. 

The second contribution H^2{'^ = ^^^ ) describes scattering terms between the 
spin waves and itinerant electrons. It reads: 

^X2 = -^Y1 ^la^P+r-q,a [{'^^XA - 1) Sln (r + f^ COS Q X 
X [cos Cp+r-q + {'^^XA - 1) Cr sin Cp+r-g] Vpq 

~^Y^ ^l^f^P+r-q,a [i'^^XA - 1) Sln (r + ^ COS Cr] X 
pqr,cr 

X [-Crsm(p+r~q + i'^SxA - 1) COsCp+r-g] Vpg 

~2§Y1 Pl^%+r-q,a [" Sin Cr + {'^^XA " 1) ^ COS Cr] X 
pqr,a 



X [cos Cp+r-g + {'^^XA - 1) fX sin Cp+r-g] Vpq 

Jk 

2N 



~^rYl f^LPp+r~q,a [" Sin Cr + {26xA " 1) (TCOsCr] X 



pqr,a- 
X [-0- sin Cp+r-g + i'25xA - 1) COS Cp+r-g] Vpg (24) 

where Vpg is a bosonic operator: 

Vpg = ala^Ch (^,) Ch i^p) + a_pagCh {^g) Sh (^p) 

+ alalgSh i^g) Ch i^p) + (aL,a_p + 6pg^ Sh {^g) Sh (^p) (25) 

By inserting the identity 1 = Spg + {1 — 6pg) in Eql24] we can distinguish between the 
contributions containing the diagonal and non-diagonal parts of Vpg. The term in H^2 
containing the diagonal part Vpp6pg, after symmetrizing with respect to ±g, can be 
further decomposed as H^j^^^iag + HAt32,diag where the boson-independent term 

HA'B2,,^aq = -^ E i^^ (2Cp)] fc - O (^) E Sh^ (^.) 

p,(7 ^ ^ g 

- ^ 5^ a [cos (2Cp)] {f3l^a:p^ + </?pJ (|) ^ ^^" (^'^) (26) 



p,(j 



contributes, together with H^^ + H^^ (EqJ251). to the electronic Hamiltonian, yielding 
modified band energies and an additional inter-subband Kondo-induced hybridization 
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of the bare AF electrons. The remaining term H^i^^ diaq contains bosons and reads: 

HXB24^a, = "^ E ^^^ (^Cp) fc - <.) (^\ J] Ch (2^,) (olo^ + al^a_^) 

+ T 5^ '"'' '°' ^^^"^^ (^^'^"^'^ + "^'^^^'^^ (ji) ^ *^^' ^^"^'^ (''^''^ + ''-^''- 



+ T S [~ '^^ ^^^^^] (''^"'^ - O ( ]^) E sh (2^.) («5«-. + 



%«-g 



p,a 
pq,CT ^ ' 9 

(27) 



Finally, the non-diagonal part of H^2 + -^62 i^- 

^ AB2,nondiag — IT \1^ ) Z-^ ^^^"^ ^^^ ~'~ ^^+''-9) \PraPp+r-q,a ~ ^ra^p+r-q,cr) 
^ ' pqr,a 

- a cos (Cr + Cp+r-g) (/5;^^ap+^_g,^ + aj^/5p+^_g^^) ]\/ g (1 - (5pg) (28) 



3.2. Effect of H^ on the electronic Hamiltonian. 

The appearance of the above-mentioned Kondo-induced hybridization terms between 
the itinerant electrons (Eqs. [22] and [2B1) suggests to perform a joint diagonalization of 
Hhand and such terms. It is convenient to define the number of AF spin waves at zero 
temperature for the Heisenberg Hamiltonian N'sw^ expressing the zero-point deviation 
for the local moments when J^ = , and the measurable amplitude of the local moment 
polarization (Sq) as: 

^s^ = (I) E Sh^ (^.) (^o) = S- ATZ (29) 

where "^q was defined in Eq. [9l 

Thus, the hybrid Hamiltonian to be diagonalized may be written: 



Ho' = E iE>p^ + «) + ^^^ E [^- (2c.)] fc - O 

per p,(T 

+ ^^^ E ^ [^°^ (2Cp)] (/5i.«p. + </5p.) (30) 



2 

p,a 

The diagonalization is realized by introducing the hybridized Fermi operators ApJ, BpJ 
through the unitary transformation: 

«Ja = ^U cos ^pa + Bl„ sin ^p^ 

PL = -^L sin ^pa + Bl^ cos ipa (31) 
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The diagonalization condition requires: 



tan (2^^ 



Jk{S^) cos {2Q 



pa J 



E^ -E- + Jk{S^) Sin {2Q 

Jk {Sq) \e 



(32) 



(^E^ - E-) ^$%2 + f/2(,.)2^ + f/|j^, (,.)^^| (^.) 



where in the second hne we have exphcitated sin (2^^) and cos (2(^p) . 
Therefore the electronic Hamiltonian in diagonal form reads: 

H^' = E K< + CO 

pa 

with hybridized energies given by (X = A,B): 



(33) 



-U 



1 



'XA 



^/3 



E^p - E- + Jk {S^) sin (2Cp) + [Jk {S^) cos (2Cp)] 



(34) 



The evaluation of the band AF order parameters in the hybridized basis, as detailed 
in the Appendix, yields: 



(s') 



AB 



^g^ (^^<) ^ E {^i^ [2 (C, + Q] «> - sin [2 (C, - Q] {nm 



p,a 



Notice that in the limit J^ — ^ , i.e. ^p — > we recover the result of Eq{T6]for the 
isolated band 



p,a 



n; 



pa 



>] 



(36) 



Conversely, for the case of a band, too weakly correlated to order antiferromagnet- 
ically by itself, i.e. for ?7 ^ , Cp ^ , we get the Kondo-induced band staggered 
moment 



];m {s^)^s = -sgn (Jk) ^ E ^^^ (K) KO + (^^>] 



(37) 



p,a 



Therefore, due to Eq. [321 in this case the band polarization is explicitly proportional 
to the effective field {— \ Jk \ {Sq)) provided by the local moments, though oppositely 
oriented. 



3.3. The longitudinal Kondo term H^ in the hybrid < ApJ,BpJ > basis. 

The longitudinal Kondo term obtained above consists of two contributions, namely 
^AB2,diagi EqJ271) and H^^f^^^nondiag (EqEHI) which have to be exphcitated in the electronic 
hybrid basis < Apa, BpJ >. Defining Zp = (p — ^p, we get: 
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JZ _ TTZ / (t) ^(t) ^ /((t) D(t)\ 

^d — ^AB2,diag (."per ) /^pcr ^^ ^pa^-^pa) 



J, 



^j:[-sin(2Z,)( 



2N 



n^ -n^ 

' ''pa ' ''pa 



)-acosi2Z,){AlB^, + BlA^,)] 



X 



p,a 



Y^ Ch (2^,) (alog + alga_ 



9 
J, 



+ SE[--(2^.)( 



' ''pa ' 'pa 



)-crcos{2Z,){AlB^^ + BlA^^)] 



X 



p,a 



J2 Sh (2^,) (ajal, + a^a. 



(3^ 



and 






uz / (t) ^(t) ^ ^(t) R(t)^ 

■^■^ AB2,no'ndiag y^pa ) r'pa ^ -"-pa ■> -^pa J 



J, 



-^J2l- ^i^ [(^^ + V.-«)] (^J<. V.-.,. - BlB^p+r-,,a) 



pqr,a 



- (T COS [{Zr + ^p+r-g)] (^L^; 



p+7 — g,(T 



i?,Up+.-gJ]^p.(i-<5p«) 



(39) 



.4-- The transverse Kondo term H^ in the hybrid < ApJ,BpJ > basis. 



3. 

The transverse Kondo term reads: 
1 



Hr = 2'^K 



/gA,o-=± j<^B,a=± 



(40) 



,(t) 



By expressing the local moments in terms of the spin wave operators a^ , and the 
electronic part in terms of the hybrid operators < AjJJ, BpJ >, we get: 

= ^\/f E 4. V.,-.4 [Ch i^,) C\-^ [p, q) + Sh {d,) C+^ [p, q)] (41) 



pqa 



+ T \/f ^ < A^+,,_.a_, [Sh (^,) C\-^ {p, q) + Ch (^,) C+^ (p, g; 

P<J(T 

+ T Vf 5Z ^4<. V.,-.4 [Ch (^,) C+^ (p, g) - Sh (.9,) Cl~^ {p, q 



+ 



Jk / 5* 



p<jcr 



p+gj-CT"--!? 



[Sh (^,) C+^ (p, g) - Ch (^,) C+;, (p, q)] 

(42) 



+ (A ^ 5) 

where the numerical coefficients Cxy ip^ l) {^i ^ = ^)B ) , depending on the angles Cp 
and ^p (EqsilS] and [32]) are explicitated in the Appendix. 
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4. Perturbative derivation of the effective magnon Hamiltonian. 

We will now describe the perturbative treatment performed to derive the effective 
second-order Hamiltonian for magnons. We begin by rearranging the total Hamiltonian, 
separating it into a basic part Hq, and a perturbation I. The basic part consists of 
the bare magnon part Hneis (Eq- flOl) plus the diagonalized electron Hamiltonian Hq 
(Eq j33!) . The perturbation includes the full transverse Kondo coupling term (J-*-) and 
the boson-dependent part of the longitudinal Kondo Hamiltonian (/^), i.e. the terms 
not included in Hq . Explicitly we have: 

Ho = Hueis + Hq 

I' =Id+IL (43) 

where /J, /^^ and I^ are respectively given by Eqsl38], 1391 and H2l 

In the following , the effect of the total perturbation I = I^ + 1^^ + 1^ will be taken 
into account through a Frohlich-type of truncated unitary transformation [25]. We 
determine the generator 7?. of the appropriate canonical transformation by eliminating 
from the transformed Hamiltonian the first order term in the perturbation /. To this 
aim, we impose / + i[7l,Ho] = 0. Introducing the notation TZ = TC^ + 7?.^^ + 7^-*-, 
we decompose this constraint into three separate equations, which can be solved [26] 
yielding: T^^^'-*-^ = \i'mt^Q^^_^l'^^'^\x)dx. Each term in the perturbation produces 
a corresponding term in the generator, namely, from /J and J^^ we obtain the 
"longitudinal" generators 7?.^ and 7^^^ while from J-*- we obtain the "transverse" 

generator 7^-*- . The terms 7?.^, 7?.^^ and TZ-^ are detailed in the Appendix. By this 
procedure. [26], we obtain the second-order effective Hamiltonian for the magnon- 
conduction electron system as: 

Heff = H' + \[n:, + 7^:, + 7^^ r, + 1^, + i^]+o (/,^ /:„ i^f (44) 

Let us anticipate here that one finds that the perturbative parameter which actually 
controls this expansion, is the ratio \Jk/ Jh\ weighed by coefficients depending on the 
electronic band structure and filling, whose expresssions are detailed in the Appendix. 
These electronic coefficients effectively reduce the magnitude of the perturbative control 
parameter from the raw value \Jk/Jh\) leading to a smooth convergence of the 
perturbative expansion even when | Jk/ Jh \ is near or exceeds one. This will become 
clear when we present our numerical results for the renormalized magnons in the next 
section. 

Finally, we make a projection onto the AF fermion wavefunction to obtain a second- 
order effective Hamiltonian for the magnons Hgjr. Let us remark that, when taking the 
average (Heff) p^j.^^oveT the AF Fermi wavefunction, we find 

\ r^d(nd)^^ \ /Fermi ^ \ l."^ ' ^d{nd)\ /Fermi ^ \\^d^ Kd\) Fermi = vJ^ndi ^d\) Fermi ('^ 
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so that the effective spin wave Hainiltonian has the simpler form 

^SW = \^eff /Fermi ^ 

fccr q 

+ 1 {IK, I',] + m„ u + [^^ i^] >^,_ (46) 

where 

(nl) = [exp {£t - /i) /knT + l] '' (47) 

Each one of the perturbative contributions above can be decomposed as a sum of 
harmonic and anharmonic terms: thus 

\ ([Kd, i:d\)Ferrm = ^ E ^ (^.'^ + ^D {<% + «-.«-.) 

+ i E ^ (^r + ^r) (4^-9 + ^g^-g) (49) 



and 



-(^7^^/i), . = iy y T'''' (ala^ + a\a 

2 \ L ' J / Fermi ^ / ^ / ^ 1 y Q Q ' -Q - 

q X,Y=A,B 

+ iE E i^r+^Di^^H. + -<.--.) (50) 

The numerical coefficients entering EqsJH] - [50] are given by long and complicated 
expressions, which we detail in the Appendix. Here we just point out that both the 
harmonic and the anharmonic parts in H^{y have contributions from both longitudinal 
and transverse Kondo terms. 

By further defining 

n$, = h^^, = ^J"'^ + i (T/^ + T/^ + T/^ + r/^) + ^ {hv'^+ + hv^) 

and 



X,Y=A,B 

we arrive at the overall effective Hamiltonian: 

(H) Fermi = E^(^ + *'?j («SS + «-g^-«) + E ^^^ (^^-9 + «g«-g) + ^0^) 
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With one last Bogolyubov transformation 



4 = 4<^h M + ^-gSh i%) 



where 



Th (2r/J 



h^r^ 



nj2 + % 
we diagonahze the effective Hamiltonian for the spin excitations (Eq j53!) . yielding: 



(54) 



(55) 



h: 



eff 
SW 



where: 



n„ 



>> 

E« 

9 






sgn(fi, + 2$,)7(fi, + 2<l>,)'-4^2 



(44 + d^-qd. 



(56) 



(57) 



is the renormalized frequency of the antiferromagnetic spin waves. In almost all the 
cases investigated numerically we have found that <l>q < so that the overall effect of 
the interaction of the local moments with the AF band is, in general, a softening of the 
spin waves with respect to the non-interacting case. Hardening for some wave vectors 
was obtained only for extremely large values of | J/f| ~ t, outside the range of Kondo 
coupling values estimated for the heavy fermion compounds addressed here. Notice 
that, in the absence of the effective anisotropy field Ba produced by the crystal field, 
at g = 0, due to limg^o^g = 0, one would get fig < 0. A non-vanishing crystal field 
thus appears necessary for the stability of the renormalized spin waves. Also, our 
result of EqjST] suggests that the gap measured at g = should not be taken for a direct 
estimation of Ba, because it also depends on the value of \Jk/Jh\- 

It is interesting to mention that the observed number of Kondo-renormalized spin 
waves in the interacting system J\fsw is: 



-) y (dM = (- 



+ 



E 

Q 

2 

N 



\^q^q 



>Ch(2ry,)+('lWsh^(ry, 
>t6j Ch (2^,) Ch (2r/,) 
Y, [Ch (2^,) Ch (2r^,) - 1] 



(5^ 



thus being always larger than the bare magnon value A/^^ (Eq l29|) . At zero temperature 
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we find that the observed Kondo-renormahzed local-moment polarization (5*^) is: 



lim (S|,) = S - i (I W |Ch (2tf,) Ch (2,,,) - 1] 



^-U|)E 



\nj2 + ^. 



V^^ J{QJ2 + $,)^ - ^l 



(59) 



indicating that the screening of the local moments due to their Kondo interaction 
with the electrons enhances the zero-point-motion quantum fluctuations of the local 
moments, already present in the bare AF Heisenberg case. This allows for a physical 
interpretation of the softening effect which, as anticipated above, we obtain for these 
renormalized magnons (detailed results in next section). The AF arrangement of local 
moments produces an effective magnon Hamiltonian in the form of a harmonic plus 
an anharmonic term, which we reduced above to a simple harmonic oscillator form by 
the final Bogolyubov transformation. That type of transformation entails the increase 
of the zero-point motion of the effective local moments, thus reducing the effective 
local magnetic moments with respect to their value in the absence of Kondo coupling. 
The relevant scale for the bare Heisenberg energies is determined by zJ^Sl+Ba). In 
the interacting case, by diagonalization we arrived to an expression of the effective 
spin-wave Hamiltonian (Eq. 1561) representing it as a new Heisenberg-like harmonic 
Hamiltonian in terms of an effective local moment < S^ >, reduced with respect to 
the full local moment S. So now the relevant energy scale for the renormalized spin 
waves is z < S^ > Jui+Ba) , a value lower than in the non-interacting case due to 
the Kondo screening of the local moments. In other words, the dynamical screening of 
the local moments is naturally reflected in a Kondo-induced softening of the magnon 
energies. 



5. Results and discussion. 

In previous section, we obtained a formally simple final expression for the renormalized 
antiferromagnetic magnons given by Eq. fl57j) . but one which depends on a series 
of coefficients which are detailed in the Appendix. The quite complicated explicit 
expressions for these perturbative coefficients, depend on the combined effect or interplay 
of the different model parameters, and most coefficients involve multiple summmations 
over the reduced Brillouin zone. Therefore, our renormalized magnon results can only be 
evaluated numerically, exploring wide ranges of the different model parameters, in order 
to assess and compare the main effect of each of them. We will show that our model can 
reasonably explain experimental magnon results [2011211 [23] employing parameters in the 
range which has been independently shown [27] to be appropriate for a phenomenological 
fit of specific heat measurements in this family of antiferromagnetic heavy fermions. Our 
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Figure 1. Simple cubic lattice: 1st Brillouin zone and special symmetry points. 



exploration of wide parameter ranges has also allowed us to verify that the convergence 
radius of the perturbative series which determines the magnon renormalization is much 
wider than one might naively have expected. As mentioned below Eq. fH4|) . the results 
presented in this section evidentiate that the actual "small parameter" controlling this 
perturbative expansion is not just the bare \Jk/Jh\ ratio. Indeed, in the renormalized 
magnon frequency this quantity appears always weighed by electronic structure- and 
filling-dependent coefficients which, when combined with suitable values of the other 
model parameters, effectively reduce the control parameter value from this ratio, leading 
to convergent results also for \Jk/Jh\ > 1, as will be shown in this section. 

The numerical study of our model has been done assuming two simple cubic 
interpenetrating magnetic sublattices, for simplicity and without much loss of generality 
(rigorously, under this assumption our results would correspond to a cubic (bcc) 
chemical lattice). Notice that CeIn3,[2T] one of the compounds where inelastic neutron 
scattering (INS) on single crystals has measured the magnons we aim to describe, is 
cubic (though fee) and has a three-dimensional Neel-type antiferromagnetic structure. 
We evaluated the magnons at zero temperature, assuming an underlying 3D Neel- 
type antiferromagnetic ground state of the system. Our study focuses on parameter 
sets far away from the quantum critical region of these systems, i.e. deep inside the 
antiferromagnetic phase. In fact, our parameters lie well inside the AF stable region 
recently determined by a DMFT + NRG study[28] of the magnetic phase diagram of the 
correlated Kondo-lattice (corresponding to the Jh = case of our model: the addition of 
non- negligible AF-like RKKY coupling Jh, as done here, will only increase the stability 
of the AF phase). As experimentally observed, [21] in this range one might expect 
Kondo-type spin fluctuations to be less relevant, and the dispersive spin waves, object 
of our study, to appear in the AF phase. For simplicity, but also in accordance with 
experimental indications [21] , we have further assumed that there is one isolated crystal 
field level of Ce^"*" which is relevant, hosting a spin S = 1/2 (in fact, S = O.S/i^ is the 
experimental magnitude of the local moments in CeIn3[2T] ) . The hopping parameter 
t was taken as unit of energy, being W = 12t the total bare electron bandwidth. 

We shall start by discussing the general trends we have found in our numerical study 
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of the renormalized magnons given by Eq. EH to subsequently focus on the description 
of the measured magnons in antiferromagnetic heavy fermion compounds. Once one 
takes into account the relative orientation of the band and local moment staggered 
polarizations according to Eq{T71 the results of our second-order perturbative treatment 
turn out to be independent of the sign of Jk- 

Fig.[T]depicts the first Brillouin Zone (BZ) of the simple cubic lattice, and we include 
the notation for the special symmetry points and BZ paths on which the spin waves 
were numerically evaluated (F = O denotes the zone center, A = X = O.Svr/a (1, 0, 0), 
Y = 0.57r/a(0, 1,0)),Z = O.Svr/a (0, 0, 1) ). Let us mention here that for the BZ 
summations we have used the special-points BZ sampling method by Chadi-Cohen (CC), 
at 4th order, for the simple cubic lattice. [29] To obtain dressed magnons with the correct 
symmetry of the lattice, and to take into account that the multiple summations involve 
linear combinations of wavectors, we have noticed [26] that it was not sufficient to 
include the basic CC set of wave vectors of the (reduced) first BZ octant, but one needs 
to extend it to the full first Reduced BZ (RBZ). Therefore one has to consider all the 
vectors obtained by applying the 48 symmetry operations of the Oh group to the basic 
CC set. Thus, at 4th-order of the CC method, we have included 5760 special symmetry 
points for each RBZ summation. To achieve higher accuracy for our determination of 
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Figure 2. Hybridized (A,B) electron bands (t) along selected BZ paths. Parameters: 
S = 1/2; T = 0;t = leV; Jn/t = 0.001; Jr/Jh = 20; n = 0.999. U/t values as detailed 
in inset. 
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the Fermi level (a delicate issue close to half-filling), we have summed over the 5th order 
symmetry-extended CC vector set, thus having used 39168 special points. 

In Fig. [2] we show the typical conduction electron bandstructure near half-filling, 
given by Eq. (1M|) . Notice that for each value of U only the lower band (denoted 
A) is filled. Chemical potential values obtained for the cases shown, are: fi/t = 
—0.63,0.84,2.30,2.33,1.13 (the Fermi level is slightly below the top of the lower 
subband), respectively for U/W = 0.008,0.25,0.48,0.49,0.67 (i.e. the U/t values of 
Fig. [21 Near the two values: U/t = 5.8, 5.8652, the conduction electrons start to develop 
an AF spin polarization (see polarization values reported in Table I), which at U/t = 8 
has increased to 0.35. The AF solution has a direct band gap, determined by the AF 
hybridized subband energies at the BZ point R on the cube diagonal (see Fig. [1]). The 
size of the gap, as well as the energy of its centroid, increases with the magnitude of 
the electron correlation U, while also a correlation-driven band-narrowing effect is seen 
to appear. For small Jk , the AF band gap value essentially depends on U (s)^^ , with 
(^)ab ~ (^)ai3 gi'owing with both the band filling n and the correlation U. In Table 1 
these trends of the band polarization with the different parameters are evidenced. 

It is also interesting to compare our band polarization values (tabulated in Table [I]) 
with those reported in Fig. 2 of Ref. |28j, for corresponding parameters. Notwithstanding 
the different respective treatments for the band electron correlations, we find quite 
reasonable agreement where we could check it. In the present work, focused on 
applications to the AF heavy fermion compounds where magnons were measured, we 
have used antiferromagnetic Kondo coupling values in a relatively narrow region around 
Jk/W ~ 10~^ — 10~^. In agreement with Peters and Pruschke, [28] iii this parameter 
range we find spin polarizations characteristic of the RKKY regime: the local moments 
are almost fully polarized, while the corresponding band polarization obtained for f/ ~ 
(i.e. unpolarized "PM" bare conduction band, corresponding e.g. to our data of the 
first two columns of Table I) is proportional to the "effective field" provided by the local 
spins, as mentioned in the last paragraph of Section 13. 2[ Our U/W = results thus 
agree quite well with the corresponding ones of Ref. [28] Since their results are given 
only for U/W = 0, 1 we could not make the comparison for intermediate U/W values. 

We have also explored wider ranges of values for Jk, to verify the consequences 
of Eql59l We have found that, if the renormalized spin waves are real and positive 
everywhere in the RBZ, then the local moment value is scarcely affected (reduced). 
Conversely, the larger are the regions of the RBZ where the EqJST] yields either negative 
or imaginary results, the stronger is the reduction of (S^), in qualitative agreement 
with the results of Ref. p8] for AF Kondo coupling. 

In Fig. [3] we show the dependence of the renormalization of the antiferromagnetic 
magnons on the bare Jk/Jh ratio at half- filling ( the most relevant filling for AF heavy 
fermion compounds) when the correlation U = 3t is not strong enough to appreciably 
AF-polarize the bare band [U {s^)a3 ~ ): we labelled this case as PM, corresponding 
to a paramagnetic bare conduction band (at Jk = 0). The top curve represents the 
bare magnons Qg (independent of the conduction electrons). We have included a small 
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Figure 3. (PM) Jk I Jh dependence of renormalized AF magnons: magnon energy 
along selected BZ paths. Parameters: S = 1/2; T = 0;i = leV]JH/t = 0.001; n = 
0.999; U /t ~ 6; Ba/t = 0.0005. Jk/ Jh values as detailed in inset; wQ denotes the bare 
magnons fig. 



anisotropy field Ba-, as inelastic neutron scattering experiments reporting a magnon gap 
indicate. [201 EI] As a general trend, anticipated in the last paragraph of previous section, 
we find that the renormalization reduces the spin wave frequency Vtq with respect to 
bare f2g, the effect growing with the ratio \JkI Jh\ ■ The softening is present throughout 
the whole RBZ, and is strongly dependent on the wavevector, being maximal around 
the r point. In Fig. H] we show the same quantities when the correlation [/ = 6t is strong 
enough to start polarizing antiferromagnetically the bare conduction band, see Table I 
(therefore we have labelled this: AF). For higher U /t one finds that the convergence of 
our perturbative series for the renormalized magnons (leading to physical non-negative 
energies) is limited to a more restricted range of values of the ratio \JkI Jh\, with respect 
to the low correlation case. The reason for this behaviour is that the fermions which 
effectively interact with the local moments are those in a neighbourhood of the Fermi 
energy of width ^ hQqp, as can be seen by looking at the explicit expressions of the 
renormalization coefficients in the Appendix. We use a Frohlich-type of transformation 
as in the BCS theory of superconductivity, so that the same type of considerations about 
the effective interactions apply. Though not too visible in the cases shown of Fig. [2l 
which lead to well-defined renormalized magnons (except U/t = 8), we have checked 
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Figure 4. (AF) Jk / Jh dependence of renormalized AF magnons: magnon energy 
along selected BZ paths. Parameters: S = 1/2; T = 0;i = leV^Jn/t = 0.001; n = 
0.999; U /t ~ 6; Ba/t = 0.0005. Jk/Jh values as detailed in inset; wO denotes the bare 
magnons flq. 



which are the main band-structure changes at larger U/t values. Near half-filling, with 
the Fermi wavector at R, for U {s^)^^ ~ {U/t = 0.1,3.0) the AF subbands disperse 
strongly around Ep resulting in a weak interaction with the local moments. But as the 
correlation increases, and U {s^)aB too, the lower AF subband, which contains the Fermi 
level around R, progressively fiattens so that now more electrons effectively interact with 
the local moments. For the cases U/t = 6,8 one also finds that another Fermi surface 
pocket appears around M. 

Fig. Odepicts the dependence of the renormalized magnons on the correlations in the 
conduction band (f/), at half-filling: the trend is a larger softening of the magnons when 
U is increased , and we show cases where U/W ranges between 0.008 (for U/t = 0.1) and 
0.49 {U/t = 5.8652). Notice that the increase of U increases the g-dependence of the 
magnon renormalization. This results from the indirect effect which U has on magnons 
(while Jk has also a direct effect, since it appears also as explicit multiplicative factor 
of the perturbatively obtained magnon corrections). U affects magnons through the 
modifications it induces in electron bandstructure, as discussed above in connection 
with Figs. [3] and H] , and through the [/-dependence of the energy denominators 
in the perturbative coefficients which determine the g-dependent renormalization of 
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magnons (details in the Appendix). The recent more refined DMFT+NRG treatment 
of correlations in an extended Kondo lattice model [28] unfortunately does not allow 
us comparison, here, as their finite U results are presented for antiferromagnetic Jk 
outside the region of interest in our problem: their correlated AF Kondo coupling 
system is studied at much too large Kondo coupling (namely, Jk = 0.5W = U) for 
the antiferromagnetic state to remain stable, the stable phase near half filling in that 
case being the Kondo insulator with all moments locally quenched. 

We exhibit the effects of doping on the magnon renormalization in Fig. El Here the 
deviation of the renormalized magnon energies from the bare magnon values increases 
with the filling: at half-filling the renormalization is largest, due to more conduction 
electrons contributing to the renormalization of magnons by coupling through Kondo 
interaction to the local moments. Doping away from half-filling we obtain a smooth 
reduction of such renormalization effects, as one would expect. Thus, both filling 
and electron correlation do increase the renormalization effects, and we have already 
mentioned that both lead to similar increases of the spin polarization of the conduction 
band (see Table I and Figs. 5-6). When entering the filling range of values at which the 
bare electron band ( Jx = 0) develops an itinerant AF polarization, the renormalization 
effects become much stronger, as shown by the very different behaviour of the cases 
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Figure 5. U dependence of renormalized AF magnons: magnon energy along selected 
BZ paths. Parameters: U /t as detailed in inset; Ba/t = 0.0005; others as in Fig. 2. 
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Figure 6. Filling (n) dependence of renormalized AF niagnons: magnon energy 
along selected BZ paths. Parameters: n values as detailed in inset; U/t = 5.865; 
Ba/t ^ 0.0005; others as in Fig. 2. 



n = 0.9 (PM band for Jk = 0) and n = 0.999 (AF band for Jr = 0). 

Let us briefly refer again to the g-dependence of the AF magnon renormahzation we 
find. Some anisotropy is present: a larger g-dependence is noticeable along BZ diagonal 
paths such as O — M or O — -R (see e.g. Fig. [6]) or paths along the symmetry axes. The 
renormahzation effects are more pronounced at long wavelenghts: in particular, they 
are maximal at the BZ center where we find a spin stiffness decreasing with doping, and 
increasing with U or Jr/Jh- Making allowance for the quite different systems involved, 
let us mention that the renormalized AF magnon behavior we obtain contrasts with the 
one recently disclosed by INS measurements in ferromagnetic metallic manaanites:^T\ 
where at low-g the spin wave stiffness appears insensitive to doping, while magnons 
exhibit a doping- dependent renormahzation at the BZ boundaries (recently suggested 
to be related to electronic correlations [32]). 

At this point, let us compare our results with the few INS magnon measurements 
available for single crystals of antiferromagnetic heavy fermions. Comparison in more 
detail may be made only with Celna, which is cubic (though f.c.c.) and presents a three- 
dimensional AF order as we have assumed for our calculation. The sets of parameters 
we find as allowing us a reasonable description of the INS magnon results which were 
described in the Introduction, [21] as evidenced by inspection of the cases presented in our 
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figures, are similar to the parameter ranges independently suggested by other authors 
for this family of compounds. A concrete example is the fit of experimental specific 
heat curves for Celns , made by Lobos et al.,[27] using a model related though not 
exactly identical to ours, within a phenomenological approach, who find: Jn/t = 0.0014, 
t = 0.5e\^, n = 1 and Jr/Jh = 980. For CeRh2Si2 they instead estimate: Jn/t = 0.0034 
and Jr/Jh = 430; and their data extrapolation for CePd2Si2 was: Jn/t = 0.0034 with 
a negligible Jr/Jh-^I\ 

6. Summary 

In the present work, we have studied spin wave excitations in heavy fermion compounds 
with antiferromagnetic long-range order , where a strong competition of RKKY and 
Kondo screening is present, as evidenced by nearly equal magnetic ordering and Kondo 
temperatures. We have described these systems using a microscopic model including 
a lattice of correlated /-electron orbitals (as in Ce-, U- compounds of this family) 
hybridized with a correlated conduction band, in the presence of competing RKKY- 
Heisenberg and Kondo magnetic couplings. 

Through a series of unitary transformations we perturbatively derived a second- 
order effective Hamiltonian describing the antiferromagnetic spin wave excitations, 
renormalized by their interaction with the conduction electrons. We have numerically 
studied the effect of the different parameters of this effective model on the magnon 
energy renormalization. Apart from the expected increase of renormalization effects for 
larger Kondo coupling, we identify another relevant ingredient. Magnon renormalization 
is also amplified by spin polarization of the conduction electrons: either if it originates 
from large correlations between the carriers or by an increase of the electron filling. We 
have been able to find appropriate sets of model parameters to describe the few existing 
measurements of magnons by inelastic neutron scattering in single crystal samples of 
antiferromagnetic heavy fermion Ce compounds (such as Celna, CePd2Si2, CeCu2). Our 
parameter sets agree with the ranges independently proposed for these materials, by 
phenomenological fits of other experiments like specific heat. Our results may provide 
information of interest for the prediction of inelastic neutron scattering experiments in 
other systems of this family, like CeRh2Si2, where there have been suggestions that the 
RKKY coupling should be stronger than the Kondo effect, [15] and the only existing NTS 
measurements are of poor quality: they were made on polycrystals|15j many years ago, 
with lower resolution. 

As outlook towards related future work, we might mention the description of the 
experimentally reported magnon damping effects, and the study of the coexistence of 
antiferromagnetism and superconductivity in the context of the present model. 
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7. Appendix. 

7.1. Evaluation of the AF order parameter in the hybrid < ApJ, BpJ > basis. 

Using tlie real-space representation we have to evaluate 
1 . 1 

Notice that EqED] implicitly assumes that (sq) is positive on the A sublattice, 
which is correct only for a FM Kondo coupling. To keep track of the correct sign for 
arbitrary sign of Jk we have to use EqJTTl yielding (^Sq) = — sgn(J/^) |(sq)|. 

For a generic site m E AU B one has the standard decomposition: 

cL = ^ E i-""*" + ^ E i.a.--'"*''"'- (61) 

^ pGRBZ ^ '■^ p&RBZ 

Substituting for c^p^^c^j^Q^ the < apa, /Spl- ^operators and recalling that on the A 
sites Rm = 2ma while on B sites Rm = (2m + 1) a one has: 



N 

peRBZ 

Rj eB c]. = 7^ 5^ [Ha - PU) cos Cp-a (4, + (3l^) sin (p] e"'^^^ 

^ peRBZ 

(62) 

Expressing the < api , /Sp]/ ^operators in terms of the< ApJ, BpJ ^ basis the real space 
Fermi operators in the hybrid basis for Ri & A are given by: 

i = ;;^ E [(4. + ^l) ^°^ (^p -Q + ^ (4<. - 4^) sin (Cp - Q] e^^^^' 

p 
and, for Rj G ;B, by: 

sl = :wYl [(4^ - ^J-) ^°^ (^^ + ^p) - ^ (4<. + ^JJ sin (Cp + ep)] et^s^ 
p 

It follows: 

^ E ^^'- = i^ E 4.4^^ {^ + ^i^ [2 (Cp - ^p)]} + ^ E ^J'^^'^ {^ - ^i^ [2 (Cp - Q]} 

leA,(T p,(j p,a- 

+ 4^ E (4.^P<. + ^J<. V) ^ cos [2 (Cp - ^p)] (64) 

p,a 

and 

-^ E ^"^- = -4^ E 4-4- {^ - si^ 12 (Cp + ^p)]} - ^ E ^J-^p- {^ + si^ [2 (Cp + Q]} 

jeBa p p 

+ i]^ E ^ ^°^ [2 (Cp + Q] {AlBp^ + B^Ap^) (65) 

p 
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The EqsJM] and ESI when substituted into EqJHDl yield: 

'AB = ^Y. i^i^ [2 (Cp + Q] AlA^^ - sin [2 {(, - Q] BIB,] 

+ i^ E ^ «^P- + ^J-^P-) {^os 12 (Cp - ^p)] + cos [2 (Cp + Q]} (66) 

When taking the Fermi average of Eq i66l the terms A^ B + B'^ A do not contribute 
and we obtain Eql35] 



7.2. The coefficients Cxy {k,q) in Eq\42- 
It is convenient to define for short 

Ll^ = cos {Z, + AZj M^^ = sin {Z, + AZj (67) 

A = ± Zp = C,p — S,p 

notice that L^^, = cos (2Z^) and M^^, = sin(2Z;,). In Eq|l2]the coefficients Cl^y {^il) 
and C^Y {k, q) with X,Y = A^B then read: 

C^xY (k, q) = 5xY [Ll,^, + {25xA - 1) M,-;,^J 

+ (1 - 5xy) [Ll,^, + (1 - 25xb) M+,^J 
C^+ (A;, g) = <5xyC+^ (A;, g) - (1 - ^xy) C+^ (A:, q) (68) 

y.c?. r/ie generators TZ^ and TZ^^ . 

In the perturbation we have several types of contributions. In /J and /^^ we find terms 
with two types of products of Bose operators. The generator corresponding to the 
perturbation term of first type (number-conserving) like YlipqT a ^pq,'^-^r,a'^p-q+r,a^p^q i 
where X = A,B while Cpqa is a numerical coefficient, is given by: 

V ^ ^pqa 



p„j. ^ ^r,a ^p-q+r,a + (i ("p "q 



ra p—q+r,crpq 
pqr. 



-^ra^p-q+r,<j^p^q l^yj 



The generator corresponding to the perturbation term of the second type like 

}-^pq,CT^Pil,cr^p,cr^p,cr^<l C'-gfeadS 



jpq 



_ q "--q 

0„. ^ O^ ^ —1— Zj/ til t/rt 

pg,CT 



•nn.rr ' '" J^'" ^ 



Xt y a^aL^J (70) 



where the (+) sign applies for bosonic creation operators, and (— ) for destruction 
operators. 
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7.4- The generator TZ^ . 

The generator TZ'^ resulting from the transverse Kondo term /^, Eql42| can be written 
as the sum of four contributions: TZ^ = ^^ y=a b '^xy ; where: 



kqcr 



)a]x 



X 



and the coefficients W^^ and 2^^ are given by: 

[Ch (^,) C+y {k, q) + Sh (^,) C3^+ (fc, q)] 

[Sh (.9,) C+^ (A;, g) + Ch (^,) C^+ (A;, g)] 



Wf/ 



-^XY 



(^f - ^r+. - ^q) 



(71) 

(72) 
(73) 



7.5. The coefficients o/ (1/2) ([7^, J]) 



Fermi ' 



We have obtained the effective Hamiltonian in Eqj56l Taking advantage of the electron- 
hole symmetry, the coefficents will be now explicitated assuming the paramagnetic band 
filling per site ra < 1 so that (^f g.) = in the ground state. 



7.5.1. The coefficients of (1/2) ([7^^, /J]) 



Fermi 



In EqHHlthe coefficients read: 



^«" = -#(1) ' E K,Kr K + (1 - «..) («4>] {<) ^^ 

^ ^ pra '' 



(74) 



+#(^)5:w,r 



per 



+ 



£A _ CB ' £A _ CB _ OhO 
L^pa '-'pa '-'per '-'pa ^m^q 



«> Sh= (2tf,) 



+ 1 (I) E (i?,) 

^ ' per 



«> 



CA _ CB 
'-'pa '-'pa 



and 



72/9 

-•anhar '^ K ' 



IQ \N 



E ^.tp^^r. [K + (1 - ^^'p) (^;t>] «> ^ 



(75) 



pra 



+#(^)E(i;,r 



pgo- 



+ 



£A _ CB ' CA _ CB _ ohQ 
L'-'pa '-'pa '-'pa '-'pa '^ib^Lq 



«> 



Sh (4^q 



7.5.^. T/ie coefficients of (1/2) ([7^:„/^,]) 
and rog^*" of Eq|l9] it is convenient to introduce: 



Fermi ' 



To write down the coefficients di"^ 



ch(^, + ^p) = c:,p sh(^g + ^p) = 6 



qp 



(76) 



Spin wave excitations in the antiferromagnetic Heisenberg-Kondo model for heavy fermions. 21 
By defining 



pqr r,p—q+r \ q,p—q+r,p ' —p,p—q+r, — qj QP \ p~q+r,(j / 

-T+ (^^^^ ^ y^B2 w /^^ \ 

■'^r,p-q+r \'^q,p-q+r,p "•" '^-p,p-q+r,-q! ^QP \"'p~q+r,(j / 

■X- A/f+ ( y^^'^ -I- y^^'i \ p^ In^ \ 

T ivlj,^p_^j_j, \^<^qp_g_i_j,,p T '^_pp_qj_j.^—q) ^^qp \"'ra / 



BA3 



■'r,p—q+r 



( yBj 



_l_ yBA-i \ (2^ /A \ 

q+r,p ' '^-p,p-q+r,-q) ^~'gp \'''rcF / 



and 



C 



pqr 



+M, 



X 



AAl 

-p,p—q+r, 



I ■yAA2 
-q ' q,p—q+r,p 



^-qp yl'rcr) 



+ L. 



r,p—q+r 

( vBAl I -yBA2 \ rt- /„A \ 

r \^-p,p—q+r,—q "•" '^q,p—q+r,p) ^qp \"'rcj / 



r,p—q+ 



M, 



r,p—q+r 



(-V, 



AA4 _|_ -yAAA 

q,p—q+r,p ' —p,p—q+r,—q 



) *^qp \1T'p-q+r,a/ 



A 
P- 
A 



'^r,p—q+r \'^q,p—q+r,p ' '^—p,p—q+r,—qj ^-^QP \'''p—q+r,a/ 



where {X = A, B 



X^I} = - f^l^^ :i^ '^" ^^^ + ^""-""A Ch (^^) Ch (^.) 



xyw 



l^y ^x+y-w ' '''V'^w "xjj 



A'ff 2 = - f^^,^^ ~y '^^ ^^' + ^---+^! sh (^^) sh (^.) 



xyw 



X 



XX3 _ (2'5xA — 1) sin {Zy + Zx_w+y, 



xyw 



Sh(^JCh(^J 






XX4 _ (25XA — 1) sin {Zy + Za;_^ + jy 

x+y—w il^\}-l'w + ": 



xyw 



-X 



X, 



XYl 

xyw 



X. 



XY2 
xyw 



X. 



XYZ 

xyw 



X. 



XY4: 

xyw 



COS [Zy + Zx-w+y) 
\fX _ CY 4- f) (Q —Q 

COS [Zy + Zx-w+y) 

\fX _ cY _|_ ft fO _ O 

I'-'y '-'x+y—w T^ "-^^^'ui ^ 'a; 

COS [Zy + Zx-w+y) 

\Cy — O^^y^y, + M (i 'tU + i 'X 

COS (z;j^ + Zx-^w+y) 
\Oy Ox-^-y—w ri[\L^ -\- \Lx 






we can write 






± 



(77) 



78) 

79) 
80) 



82) 
83) 
84) 
85) 
86) 

87) 
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Next, by defining 

J ipqy ^'-^r,p—q+r \ q,p—q+r,p ' —p,p~q+r, — q) ^^qp \ p—q+r,cr / 

— /'+ (^^^^ A_ y^B2 \ p^ I A \ 

^r,p~q+r \'^q,p—q+r,p "•" '^-p,p—q+r,—q) ^^QP \'''p~q+r,a / 

' ''■^-''r,p—q+r y'^q^p—q+r^p ' '^—p,p—q+r, — qj ^QP \"'r(T / 

"•" ■'^r,p~q+r \'^q,p—q+r,p ' '^—p,p—q+r,—qj ^qp \"'tcj/ 



and 



•'^'■pqr '^'^■^r,p-q+r \'^-p,p-q+r.-q ' '^q;p-q+r,p) ^-'qp \''-rcr/ 

-U /'+ (y^^'^ -i_ yBA2 \ p^ /n^\ 

' ^r,p-q+r \'^-p,p-q+r,-q ' '^q,p-q+r,p) ^^qp \'''ra/ 

— M+ ( Y^^'^ -X- y^^'^ \ (T Ir,^ \ 
■'■^■'■r^p—q+r \'^q,p—q+r,p ' '^—p,p—q+r, — q) ^qp \"'p—q+r,(T / 

- T+ (T^ABi i yABA W /„A \ 
■^T,p-q+r \'^q,p-q+r,p "•" ^-p,p-q+r,-q) ^qp \"'p-q+r,a / 



we can write 



^<=f (I)'e(i-^-)-^ 



± 

pqr 
pr 



(89) 
(90) 



7.5.3. The coefficients of (1/2) { [7^^, I^] )^^^^. . In Eqj50], if n < 1 , so tliat (n^^) = 
in tlie ground state, one finds T^^^ = S^^ = 0. The other X"^^ coefficients, by defining 

{fk,k±q,a) = {^i-a) " (^fcig.a) , ^Cad: 

ka 
ka 

+ ^I1 ^^-.Ch (^,) C+:^ (k - q, q) (/,^,_,,,> 

ka 

fc(T 

^/^ = -^ E <f Sh (^,) C+;, (A: + g, -q) (n^.) 
+ ^ E >^M Ch {^q) C^s (k + q, -q) (nt^) 

k<T 

- ^ E ^^'l^.Ch (^,) C+;^ (fc - g, q) «_^> 

fccr 

+ ^ E ^^^^.Sh (^,) C+^ (A: - g, q) «_> (92) 
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and 

k 

- ^ E Kt^^ (^'^) ^BA (k + g, -g) (^t,,.> 

+ ^ E ^^"^-.Ch (^,) C+^ (A; - g, g) (n^,,,) 
fc 



2iV 



The 5^^^ coefficients read: 



'^.^^^ = +^ E ^A:!-.Ch (^.) CX], (k - q, q) (iV,,_,,.) 



2N 

k 



+ W E ^fcf .Sh {§,) C^B {k - g, g) (iV,,_,,.) 



2iV 



?ABl _ -^K^ S^^^AiAB , .^ ,„^ ,.. . ,,„ ^^,,^. 



'^r ' = -^ E >^fc .Ch (^,) Cl-^ {k - g, g) 

k 

+ ^ E <^.Sh (^,) Ci^ (A: - g, g) «_,) 



^r ^ = +^ E ^^^'^^.Ch (^,) C+^ (A: - g, g) (nt,,) 

k 

- ^ E ^^^'^".Sh (^,) C+;^ (A; - g, g) {nt,,) 

k 



(93) 



(94) 



(95) 



(96) 



The coefficients 5^^^ can be obtained from Sf^^ by interchanging W^^ and 



Ch ['dq) respectively with Z^_ and Sh {d 
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Table 1. Conduction band AF spin polarization, < s^ >ab, dependence on other 
model parameters, for the cases depicted in respective Figs. 3 to 6. In first column 
we enter the relative Kondo coupling magnitude (cases of Fig. 3: notice that 
Jk/W ^ Jk/Jh X 10~^, here), while in second column we state the corresponding 
band polarization < s^ >ab obtained; similarly, for Fig. 4 and the related 3rd and 4th 
columns of this table. In 5th column we enter the electron correlation U /t cases from 
Fig. 5, and next column states the corresponding < s^ >ab- In 7th column we enter 
filling n values of Fig. 6, and in last column the corresponding < s^ >ab values. 

Jk/Jh (PM) < S' >AB Jk/Jh (AF) < S^ >AB U/t < S' >AB U < S^ > AB 



0.1 


-6.10E-7 


1.0 


-7.20E-2 


0.100 


-1.14E-4 


0.100 


-3.81E-6 


1.0 


-6.10E-6 


4.0 


-7.20E-2 


1.000 


-1.15E-4 


0.500 


-2.96E-5 


5.0 


-3.0E-5 


5.5 


-7.20E-2 


3.000 


-1.22E-4 


0.600 


-3.97E-5 


10.0 


-6.10E-5 


7.0 


-7.20E-2 


5.000 


-1.42E-4 


0.700 


-5.24E-5 


15.0 


-9.15E-5 






5.500 


-1.51E-4 


0.800 


-6.92E-5 


20.0 


-1.22E-4 






5.800 


-3.52E-3 


0.850 


-8.00E-5 


25.0 


-1.50E-4 






5.8650 


-8.69E-3 


0.900 


-9.94E-5 


30.0 


-1.80E-4 






5.8652 


-8.71E-3 


0.999 


-8.69E-3 



